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Abstract: Linear rate equations are used to describe the cascading decay of an ini- 
tial heavy cluster into fragments. We consider moments of arbitrary orders of the mass 
multiplicity spectrum and derive scaling properties pertaining to their time evolution. We 
suggest that the mass weighted multiplicity is a suitable observable for the discovery of 

Q_i' scaling. Numerical tests validate such properties, even for moderate values of the initial 
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mass (nuclei, percolation clusters, etc.). Finite size effects can be simply parametrized. 
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In this work we consider binary fragmentation processes where any fragment with 
mass number k breaks into fragments with mass numbers j and k — j, j = 1, 2...k — 1, 
with a probability Wjk per unit of time. It is assumed that Wjk is time independent. By 
definition, Wjk = if j > k and Wjk is symmetric if j is replaced by k — j, naturally. For 
technical reasons, we will use for Wjk the double of the actual transition rate whenever the 
special case Wj^j occurs. Let Nj(t) be the multiplicity of fragment j at time t in a process 
initiated from the decay of a cluster A, namely iVj(O) = Sja- The model under study is 
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described by the following set of linear, first order differential equations, 
m >' ' A 

X ; ^ = -cjNj + J2 w jkN k , j = 1, -A, (1) 

d ; k=j+i 

with 

i-i 

c 



With components Nj, j = 1, ...A, for a column vector |jV >, the system, Eqs.(l), boils 
down to d\M > /dt = W\Af > with a triangular matrix W. The general solution of Eqs.(l) 
is obviously a sum of exponentials whose rates of decay in time are the trivial eigenvalues 
of the triangular W, namely the diagonal matrix elements — Cj. 

This matrix W has a remarquable property, namely a fixed left (row-like) eigenstate 
A^i, whose components M.\j = j, j = 1, ---A, do not depend on the wWs. This comes from 
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the symmetries of W demanded by the conservation of the total mass Mi = J2j=iJNj, 
with dMi/dt =< Mi\W\Af > . The corresponding eigenvalue is, naturally, —c\ = 0. 
It is thus convenient to define the "mass weighted multiplicity" (MWM) vector U with 
components Ui = iN{, whose evolution is governed by a matrix V, with matrix elements 
Vjfc = jWjk/k, hence 

dU/dt = VU. (3) 



The surface under the histogram defined by U is therefore time invariant, X)i=i Ui/A = 1. 
This invariant normalization gives a hint that an analysis of U might be the best way to 
reveal scaling and even universality properties of the fragmentation process. 



For that analysis, we first consider a moment of arbitrary order q 

A A 
3=1 3=1 

where the exponent q is any real number. The time derivative of such a moment is 
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dM q 
~df 
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E-H ^2 WjkNk-cjNj] = Y,k q - 1 U k d(q,k), 

3=1 \k=j+l I k=l 



(5) 



with 



fc-i 

d(q,k) = ^w jk 

3 = 1 



(6) 



where we have used an interchange of indices j and k in the double summation. 



The structure of the right-hand-side of Eq.(5) gives a special interest to those cases 
where d(q, k) shows scaling properties with respect to k. For instance 1 ) if a condition 
d(q, k) ~ d(q) is met, this approximate independence with respect to k induces an eigen- 
value problem where the corresponding (eigen) moments M q are decoupled and decay ex- 
ponentially with respect to time. A sufficient condition for this is a behavior Wjk — 
k~ 1( fi{j /k), whence, when k is large enough, 
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d(q,k) = -J2^U/k) 



3 = 1 
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- ' dx <p{x) (x q - 1/2), 
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where x = j/k. This indicates indeed a limit where d(q, k) becomes d(q), independent of 
k. Accordingly, 

^ - d(q)M q . (8) 

Such a decoupling of eigenmodes is also found in the QCD theory of energy-momentum 
fragmentation 2 ^, where any moments with a (continuous) order q real and positive is an 
eigenmode, in contrast however with the present discrete representation of fragmentation 
where only special values of q can be selected 1 ). 

In this letter, we study a more general case where d(q, k) ~ k a d(q) and a is any real 
number a priori. This occurs in particular if Wjk — k a ~ x (p{j /k). In such a case the resulting 
evolution equation for moments reads 

^=J2 kq ' 1U * d ^ k ) - d («) (9) 
k=l 

Here d(q) is still taken from the right-hand side of Eq.(7). This property, Eq.(9), occurs 
in continuous systems 3 ). Our aim is to rather investigate discrete fragmentations as well. 

Let us assume temporarily that Eq.(9) is an exact result. This leads to now well- 
known 3 ) solutions with scaling properties, 

M q (t) = [s(t)] q - 1 m q , s(t) = {[s(0)]- a + auot)~ l/a , m q+a = m q (1 - q)w/d{q), (10) 

where s represents an average for the cluster sizes at time t. The recursion relation between 
the coefficients m q relates the various moments to one another via the function d(q), see 
Eqs.(7,9,10), which contains the dynamical information about the system. Finally 10 is an 
integration constant. Such scaling solutions for continuous systems are valid only if a is a 
strictly positive number. 

In our present study, the physical quantities under observation are the U^s, governed 
by Eq.(3). The translation of Eq.(10) in terms of Ui clearly reads 



tf 4 (f)«[a(f)]~7[*M*)]. 
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(11) 



where / can be deduced from d(q), but needs not be explicit in a numerical verification 
of scaling. Indeed, if one solves Eq.(3) numerically, it is sufficient to observe that, for 
different times t, the various Ui(t) show similar shapes when referred to the new scaling 
variable y = i/s(t). Naturally, for each time t, one needs a knowledge of s(t), by numerical 
derivation if necessary. 

In this letter, we report numerical results for a two-parameter class of models where 

wjk = k a - 1 [{j/k)- h + (l-j/k)- b }/2. (12) 

Such an ansatz for Wjk was already investigated with interesting results for the possible 
existence of eigenmoments 1 ). We first calculate Eq.(3) for many values of A,a,b,t, then, 
as a rule of thumb, discard both U\ and Ua- This is because an obvious peak in the 
[/-histogram exists at Ua for small fs and at U± for large fs. We are rather interested 
in identifying phenomenologically s(t) as the position of the intermediate extremum (if 
it exists and if it is unique) of Ui(t) as a function of i. The reason for this procedure is 
the fact that, as already stated, the [/-histogram has a constant surface. An intermediate 
bump, or conversely a dip, between the "source peak" Ua and the "sink peak" U\ is thus 
likely. 

Indeed, as shown by Fig. la, where A = 24, a = 0.5, b = —1, and t = 1.2, 2.2, 3.2, 
there is an intermediate bump, which smoothly moves towards lighter fragments when 
time increases. Simultaneously, the source peak at A = 24 decays, while the height of 
the fragment bump grows. The case shown by Fig. lb is rather different. There is a flat 
dip rather than a maximum, when A = 24, a = 0.5, 6=1, and t = 0.28, 0.36. The 
flip from a source peak at shorter times into a sink peak at larger times, see t = 0.44, 
is striking. Finally, for Fig.lc, where A = 24, a = —0.2, b = —1, and t = 9, 10, (and 
t = 11, for comparison) we observe three maxima, namely a source peak, a sink peak and 
an intermediate flat maximum. It will be noticed that the time scales which govern these 
evolutions seem to depend strongly on a and b. 

We now turn to the search for scaling properties. For those histograms at time t 
which show an intermediate maximum U max (t) at some mass i ma x(t), see Fig. la and lc, 
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we set s(t) = imax(t). Then the mass scale is defined as y = i/i m ax{t), and we consider 
Log[Ui(t) /U max (t)] as a function of y. For those histograms which show an intermediate 
minimum, we do the same, except that the reference for scale and normalization is taken 
with respect to the minimum, naturally. 

This allows a comparison of histograms at different times: all histograms reduce to 
when y = 1. Finally we discard from our analysis those histograms for which U{ is 
monotonically increasing or decreasing between i = 2 and i = A — 1. A few results are 
shown on Figs. 2. Fig. 2a exhibits a remarkable match of such "renormalized" histograms 
for A = 24, a = 0.5, b = — 1, and t = 0.8, 1.4, 2.0. It must be stressed that the data for 
t = 0.8 cover only the interval 0.1 < y < 1.3, while those for t = 1.4 cover only the interval 
0.2 < y < 3, and finally those for t = 2 cover the interval 0.3 < y < 4.5. This is because the 
position imax decreases considerably as a function of time. Hence the mutual continuation 
of the three curves from one another is all the more striking. Another remarkable match 
is seen on Fig.2b, where A = 24, a = 0.5, b = 0.5, and t = 0.4, 0.5, 0.6, 0.7. 

We have verified that this result does not depend on A, but only the parameters a and 
b, which must be in a range a > and b < 1. Indeed, as shown for instance by Fig. 2c, with 
A = 24, a = —0.2, b = —1, and t = 8.7, 9.0, 9.3, the overlap of the renormalized histograms 
is less convincing. Indeed, it is seriously violated in the region of small fragments y < 0.2. 
Finally, the results obtained when there is a minimum rather than a maximum give no 
evidence for universality of shapes in such cases, and strongly point to the contrary, see 
Fig.2d with A = 24, a = 0.5, b = 1, and t = 0.24, 0.28, 0.32. 

For comparison and contrast with Fig. 2a, we show in Fig. 3 the same attempt for 
mutual continuation with A = 24, a = 0.5, b = —1, but t = 1.8, 2.2, 2.6. Some success is 
obtained for y < 1, but the attempt fails for y > 1. It can be concluded that a "universal" 
shape of the renormalized histograms is likely for intermediate times, but not for longer 
times. 

It is now interesting to confront our numerical findings with the solutions in the 
continuous limit 3 ) where Eq.(9) is an exact property of the system. 
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i) scaling case: a > 0, b < 1 

In the continuous limit, one obtains the scaling behavior described by Eqs.(10,ll). If 
we assume for / a simple factorized form and perform the same rescaling transformation 
as for the distributions we get 

z(y) = (1 - b')/a' (l - / + a' logy) , (13) 

where the parameters a', b' are to be compared with those a, b of the continuous limit and z 
is the scaling function emerging from the time-dependent distributions Log[L/j(f) /U max (t)]. 
A comparison between z and Log[L/j(t) /U max (t)], see Fig. (4), shows that the rescaled 
histogram Log[t/"i(t)/t/" max (t)] is well reproduced by z, Eq.(13), with b' = 6, but a' ~ 
0.35 < a = 0.5. Indeed, in the continuous limit, the scaling function z(y),y = i/s(t), see 
Eq.(ll), is known to behave as z oc (1 — b) logy for small y and z oc — y a for large y. It 
seems that the finite size effects do not affect the scaling properties that we have thus 
observed in the considered cases, albeit possibly modifying some of the scaling indices. 

ii) shattering case: a < 0, b < 1 

In this case, scaling is approximately valid in the region y > 1 and strongly violated 
in the small y region. As a matter of fact, from fitting Fig. 2b, we find empirically that the 
parametrization, Eq.(13), gives a satisfactory agreement for a' ~ 0.1 in the region y > 1 
but fails for y < 1. Note again that a' ^ a. In the continuous limit, it is known that a 
shattering transition 4 ) takes place, whereby a part of the system is quickly transformed 
into a powder of infinitesimal constituents. In the present discrete case, the behavior of 
U for small values of y can be traced back to a similar formation of a sector of lightest 
fragments of unit mass. 

iii) non-scaling (evaporation) case: b > 1 

The absence of a scaling behavior is likely due to the presence of a minimum instead of 
a maximum in the function z(y). It is to be noted that the marginal value 6=1 corresponds 
to the field-theoretical case mentionned in reference 2 ) . However, it is not excluded that a 
partial scaling could be restored by some different empirical definition of the rescaling 

6 



procedure. Nevertheless, the fact that the values b > 1 give a strong enhancement to the 
contributions of the very small (and thus very large by symmetry) fragments at the vertex, 
see Eq.(7), is a hint towards a difference between the continuous and the discontinuous 
cases. 

Discussion and Conclusion: 

We have discovered that, at those intermediate times where a sufficient amount of in- 
termediary mass fragments has been generated (and evaporation and/or shattering are not 
dominant), the multiplicity histogram has a stable shape. This shape can be parametrized 
by simple functions and simple parameters and evolves in time under simple scaling rules. 
This result has been observed numerically for many initial masses, even moderate ones, 
and a fairly general class of models for the fragmentation vertices. The result is also true 
for several excursions out of this class, which cannot be reported here in detail. It is thus 
reasonable to conjecture that a central limit theorem, due to the semi-group, iterative na- 
ture of the binary fragmentation cascade, is at work when the intermediary mass fragments 
build a bump in the mass weighted multiplicity histogram. 

We have found that our phenomenology of discrete fragmentation shows both similar- 
ities and differences with the description 3 ) proposed for continuous systems. For example, 
similarities are found when the initial mass is very large, naturally. Differences are found 
for moderate initial masses. For instance we find no obvious scaling if the histogram has 
a minimum (evaporation), and conversely, we find a residual scaling for part of the his- 
togram in the presence of shattering. In those cases where there is an agreement between 
our description and that of the continuous limit, we find nonetheless finite size corrections 
in the behavior of the scaling curve. 

Besides the conjectured existence of a suitable central limit theorem, this work opens 
several lines of investigation. One one hand, those multifragmentation data which con- 
tain sizable numbers of intermediate mass fragments should likely be analyzed in terms 
of scaling. A (nuclear, atomic) clock is not necessary for dating the corresponding his- 
tograms, since our scheme is actually time independent. On the other hand, one should 
investigate whether binary 5 ) cascades make a sufficient model, or whether fragmentation 
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is very sensitive to many-body effects in the medium 6 ) . Comparison with the properties of 
percolation 7 ) is also in order. 
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Figure Captions 
Fig. la 

Mass weighted multiplicity Ui as a function of fragment mass i for initial mass A = 24, 
and vertex parameters a = 0.5 and b = —1. See Eqs.(3,12). Notice how the "source peak" 
U24 decreases when time t increases from t = 1.2 to t = 3.2. Simultaneously, a population 
of intermediate fragments moves into lighter ones. 

Fig.lb 

Same as Fig. la, but b = 1. In contrast with Fig. la, the source peak feeds directly the pop- 
ulation of lightest fragments and even disappears for t = 0.44. The process is reminiscent 
of evaporation. 
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Fig.lc 

Same as Fig. la, but a = —0.2. Both source and sink peaks are present, together with an 
intermediate, temporary maximum, see t = 9, 10. A shattering phenomenon is likely in 
this case. 

Fig.2a 

Scaling comparison of mass weighted multiplicities at various intermediate times. For each 
time we plot Log(Ui/U max ) as a function of y = i/i max - The parameters are A = 24, a = 
0.5, b = —1, t = 0.8, 1.4,2. Notice how a unique scaling curve emerges from overlapping 
segments of curves. 

Fig.2b 

Same as Fig. 2a, but b = 0.5, t = 0.4, 0.5, 0.6, 0.7. Notice how intermediate and lower mass 
fragments follow scaling. The population of the highest masses deviates from the scaling 
curve, but its decay nonetheless contributes to it. 

Fig.2c 

Same as Fig. 2a, but a = —0.2. Scaling is likely preserved for y > 0.8, and strongly violated 
for y < 0.3. This indicates shattering. 

Fig.2d 

Same as Fig. 2a, but 6 = 1. Absence of scaling when referred to the minimum of Ui. 
Fig.3 

Same as Fig. 2a, but larger time intervals. Evidence of scaling violation when y > 1. 
Fig.4 

Same as Fig. 2a, but comparison with the theoretical (dashed line) estimate z(y), see 
Eq.(13). 
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